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^"^ ■ Abstract 

\q ' We present high-resolution (R = 20000) spectroscopy of Hjj" absorption toward the luminous Galactic 

center sources GCS 3-2 and GC IRS 3. With the efficient wavelength coverage afforded by Subaru IRCS, 
six absorption lines of H^T have been detected in each source from 3.5 to 4.0 ^m, three of which arc new. 
In particular the 3.543 fim absorption line of the i?(3, 3) z transition arising from the metastable (J, K) 
= (3, 3) state has been tentatively detected for the first time in the interstellar medium, where previous 
observations of H^" had been limited to absorption lines from the lowest levels: (J, K) = (1, 0) of ortho-H^ 
. and (1, 1) of para-H^. 

The Hjj~ absorption toward the Galactic center takes place in dense and diffuse clouds along the line of 
£SJ \ sight as well as the molecular complex close to the Galactic nucleus. At least four kinematic components 

are found in the Hj absorption lines. We suggest identifications of the velocity components with those 
of H I, CO, and H2CO previously reported from radio and infrared observations. H^ components with 
velocities that match those of weak and sharp CO and H2CO lines are attributed to diffuse clouds. Our 
observation has revealed a striking difference between the absorption profiles of H^" and CO, demonstrating 
that the spectroscopy of H^" provides information complementary to that obtained from CO spectroscopy. 



that the spectroscopy of ^ 

The tentative detection of the i?(3, 3)' line and the non-detection of spectral lines from other J > I 
levels provide observational evidence for the metastability of the (3, 3) level, which is theoretically expected. 
This suggests that other metastable J = K levels with higher J may also be populated. 
Key words: ISM: clouds — ISM: lines and bands — ISM: molecules — Galaxy: center 
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■ 1. Introduction A search for interstellar HJ was initiated by Oka (1981) 

more than two decades ago, and culminated in the discov- 

The crucial role which the H^~ molecular ion plays in in- ery of H^ in the interstellar medium toward AFGL 2136 

terstellar chemistry was first addressed by Watson (1973) and W33A, young stellar objects deeply embedded in 

and Herbst and Klemperer (1973). H. 4 " is produced by molecular clouds (Geballe, Oka 1996). Since then the 

cosmic-ray ionization of H 2 to H^ followed by the efficient study of interstellar H3 has progressed rapidly by a se- 

Langevin reaction, ries of successful observations in dense clouds (McCall et 

jj+ _|_ jj^ _^ jj+ _|_ jj q\ al. 1999) and in diffuse clouds toward the highly reddened 

star Cygnus OB2 12, the Galactic center, and many other 
It works as a universal protonator (acid) in the efficient s i g htlines (McCall et al. 1998; Geballe ct al. 1999; McCall 
proton hop reaction et a \ 2 002). These observations have not only demon- 
ic -f X — * H9 + XH + (2) strated the ubiquity and abundance of this molecular 

ion, which was anticipated in general in molecular hy- 

for most molecules or atoms X (He, Ne, N, and 2 are drogen dominated plasmas (Martin et aL 1961)j but also 

a few exceptions). Subsequent general reactions with a reve&led an interesting enigma related to the H + chem _ 

molecule Y, igtry of the diffuge interstellar medium. While the col- 

XH + + Y — ► XY + + H, (3) umn densities measured in dense clouds agreed well with 

n 1 r , 1 1 1 1 1 • 1 j 1 those predicted by the Hj" production and destruction 

lead further to chemical networks which produce complex . , 6 51 , 

, , . , 1 ii -1 i mechanisms ol Herbst and Klemperer (1973 Geballe. 

molecules m dense interstellar clouds. ^, H „ , , _ „ , H „. f. v ' v , . 



Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan. 
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diffuse interstellar medium have been found to be orders 
of magnitude higher than estimated using canonical val- 
ues of plasma chemical constants (McCall et al. 1998a, b; 
Geballe et al. 1999, McCall ct al. 2002). In order to have 
a better understanding of this problem, we undertook an 
absorption-line survey of this molecular ion toward the 
Galactic center. This study also provides unique informa- 
tion on the medium toward the Galactic center because 
of the special characteristics of as an astrophysical 
probe. 

The absorption-line survey was designed to cover the 
maximum number of absorption lines in clouds with 
various physical conditions. The Galactic center sources 
are ideal for this purpose, since they suffer heavy visual 
extinction of Ay = 25-40 (Cotera et al. 2000), the high- 
est in the Galaxy among those obscured mainly by dif- 
fuse clouds. The line of sight to the Galactic center cuts 
through dense and diffuse clouds in the intervening spiral 
arms as well as the molecular complex close to the nucleus. 
Thus, in different physical conditions is efficiently sam- 
pled in one observation. The intervening clouds have 
different radial velocities, and the nature of the clouds 
were studied in previous spectroscopic observations us- 
ing molecular and atomic tracers from near-infrared to 
radio wavelengths. Two infrared sources, GCS 3-2 in the 
Quintuplet cluster (Nagata et al. 1990; Okuda et al. 1990) 
and GC IRS 3 (Becklin, Neugebauer 1975; Becklin et 
al. 1978) near Sgr A* have been selected, since both ap- 
pear to have intrinsically featureless spectra due to dust 
emission and are sufficiently luminous to provide good 
continuum fluxes for L-band absorption spectroscopy. 

Wide wavelength coverage is essential for a line survey. 
The observation of interstellar Hjj~ has so far been almost 
entirely limited to the three absorption lines, [R(l, 0), 
R(l, l) z , and R(l, 1)"], all arising from the lowest J = 1 
rotational levels, i.e. (J, K) = (1, 0) of ortho-Hjj" and (1, 
1) of para-H^" (note that the J = level is forbidden by 
the Pauli principle; McCall 2001). This is because high 
spectral resolution and wide wavelength coverage are of- 
ten incompatible. This limitation has now been overcome 
by a new generation of cross-dispersed infrared spectro- 
graphs, such as IRCS on the Subaru Telescope. We can 
now simultaneously observe at high spectral resolution not 
only most of the transitions starting from the J = 1 levels, 
including Q(l, 0) and Q(l, 1), but also other transitions 
starting from higher J levels. This wide coverage has led 
us to a tentative detection of the i?(3, 3) z line starting 
from the metastable (J, K) = (3, 3) level. 

2. Observations 

These spectroscopic observations were made on 
2001 June 16 (UT) using the Infrared Camera and 
Spectrograph (IRCS; Tokunaga et al. 1998; Kobayashi 
et al. 2000) with the 8.2 m Subaru Telescope on Mauna 
Kea. Subaru IRCS is equipped with an echcllc and a cross- 
dispersing grating to offer high resolution spectroscopy 
without sacrificing wavelength coverage. Figure 1 shows 
the wavelength coverage of the IRCS echelle mode in the 



3 fim region. The entire L-band (2.84-4.18 (im) can be 
covered with six standard settings, with 70% coverage in 
three settings. Thus, Subaru IRCS is a good match for 
a survey of lines scattered in a single band. The grating 
setting can be modified continuously for the best coverage 
of target lines in order to maximize observing efficiency. 

In the current experiment we employed two grating set- 
tings optimized for coverage of the V2 <— Hjj" lines in the 
3 fiva. region illustrated in figure 1. At the time of the ob- 
servations, there was high dark current in the upper two 
quadrants of the lk x lk InSb detector array, so we cov- 
ered the relevant wavelength ranges using only the lower 
half of the array. A 0"15 x 4"5 slit was used to achieve 
R = 20000 spectral resolution. The spectra were recorded 
by nodding the telescope between two points separated 
by 2'.'2 along the slit to subtract the sky emission and 
dark current. Spectroscopic standard stars were observed 
through airmasses similar to those of the objects in order 
to cancel out the atmospheric transmission efficiency. The 
details of the observations are summarized in table 2. The 
seeing was not better than 0"7 at L', and varied during 
the observing period. Spectroscopic flat frames were ob- 
tained at the end of the night with a halogen lamp in the 
telescope calibration unit installed in front of the instru- 
ment window. 

3. Data Reduction 

The observed spectrograph images were stacked and av- 
eraged without any pixel registration since both the ac- 
cumulation of the instrument flexure and the telescope 
tracking error during automated guiding was negligible 
over the typical tracking period. The average frame was 
flat-fielded with a dark-subtracted halogen lamp frame. 
Bad pixel masks were created by collecting pixels with 
low, high, or varying responses, or with a high dark cur- 
rent, based on the statistics of the halogen lamp and the 
dark-current frames. Those outlier pixels were filtered out 
before spectral extraction. 

One dimensional spectra were obtained using the 
IRAF 1 aperture-extraction package. Spectroscopy in the 
3 /jm region suffers severely from interference by the ab- 
sorption lines of molecules in the Earth's atmosphere. 
Poor cancellation of the atmospheric lines by way of di- 
viding by the spectra of the standard stars is often the 
primary source of systematic errors, which easily over- 
whelm statistical noise. We carefully examined the lo- 
cations of the Hjj~ lines listed in table 1 using a custom 
written IDL code. The code handles (1) linear registra- 
tion of the wavelength offset between the object and the 
spectroscopic standard, (2) rescaling the normalized stan- 
dard star spectrum according to Beer's law to minimize 
the airmass mismatch that sometimes exceeded 10% in 
our observation, and (3) correction of the spectral reso- 
lution by deconvolving the spectra. Wavelength calibra- 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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tion was performed by maximizing the cross correlation 
between the object spectrum and the model atmospheric 
spectrum calculated by ATRAN code (Lord 1992). 

4. Results 

The result of the line survey is shown in figure 2 and 
figure 3 for the positive and negative detections, respec- 
tively. The model atmospheric transmission curves con- 
volved to the same spectral resolution are shown in order 
to discriminate between the real spectral features and the 
residual features of the strong telluric absorption lines. 
All of the absorption lines in our coverage starting from 
the lowest J = 1 levels of ortho- and para-H+ [R(l, 0), 
Q(l, 0), R(l, l) 1 , R(l, l) u , and Q(l, 1)] were successfully 
detected, while those from J > 1 levels were all negative, 
except for R(3, 3)'. 

The Hjj~ absorption lines from the J = 1 levels show 
several discrete kinematic components. The absorption 
profiles are deconvolved into Gaussian components cen- 
tered at LSR velocities of -110 140 km s _1 , -50 - 

—60 km s _1 , and km s _1 (figure. 4). The deconvolution 
parameters are summarized in table 3. GC IRS 3 shows 
wide absorption at —20 km s _1 and a wing at +50 km s _1 
on a broad pedestal absorption. A weak absorption is seen 
at -170 km s" 1 in GCS 3-2. Note that the velocity com- 
ponents are less distinguishable at 3.668 /im, previously 
observed toward the Galactic center (Geballe et al. 1999). 
This is because the absorption feature consists of a close 
doublet of R(l, 0) and R(l, l) u transitions with a sepa- 
ration of Av = 35 km s _1 . 

The column density of H^" was derived from the ob- 
served equivalent width (W\) using the equation W\ = 
(87r 3 A/3/ic)7V|/i| 2 , where N is the column density of H^ 
in the lower state of the transition. The values of |/x| 2 
listed in table 1 have been used. The velocity-integrated 
column densities of individual transitions are summarized 
in table 1. For the absorption lines with quality to al- 
low Gaussian deconvolution, the column density of each 
velocity component was calculated and given in table 3. 

5. Discussion 

5.1. Location of Absorbing Clouds 

We discuss below possible identifications of the absorb- 
ing clouds, referring to previous observations at radio and 
infrared wavelengths s(figure 5). We will base the follow- 
ing velocity discussion on GC IRS 3 spectral lines unless 
otherwise stated, since the line profiles are sharper than 
those of GCS 3-2. 

5.1.1. —140 km s _1 to —110 km Component 

The —140 km s _1 component is attributed to clouds 
in the "expanding molecular ring" (Kaifu et al. 1972; 
Scoville 1972). The "expanding molecular ring" is a chain 
of molecular clouds orbiting around the nucleus at 200 pc 
from the Galactic center, and gradually receding from it. 
The —140 km s _1 component of GC IRS 3 appears at 
— 110 km s _1 in GCS 3-2, which is reasonable because 
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the "expanding molecular ring" appears with a less neg- 
ative velocity at positive Galactic longitude. The addi- 
tional possible peak at —170 km s _1 seen towards GCS 3-2 
could be the high velocity component discussed by Giisten 
and Downes (1981), Liszt and Burton (1993), and Yusef- 
Zadeh, Lasenby, and Marshall (1993). 

5.1.2. — 60 km s _1 Component 

A clear match of the peak velocities at —60 km s _1 
is noted between the H^" lines and the H2CO line ob- 
served in absorption toward Sgr A (Snyder et al. 1969; 
Giisten, Downes 1981), and the H 1 absorption line (Liszt 
et al. 1985). The cloud velocity is consistent with the ra- 
dial velocity of the "3 kpc arm" . The "nuclear disk" rotat- 
ing at 300 pc from the Galactic center (Sanders, Wrixon 
1973) also shows a similar velocity at the longitude of 
GC IRS 3 (Whiteoak, Gardner 1979). However, we infer 
that the contribution of the "3 kpc arm" clouds is domi- 
nant, since the component is also seen in GCS 3-2 at the 
same velocity. This interpretation is also supported by 
the absence of a positive velocity component which is ex- 
pected from the "nuclear disk" for the positive Galactic 
longitude of GCS 3-2. 

5.1.3. km -1 Component 

The km s _1 component is also apparent in the H2CO 
and H I absorption spectra of Sgr A in figure 5. The 
km s _1 absorbers are usually attributed to the "lo- 
cal clouds" within a few kpc of the solar neighborhood, 
which is most evident in the H 1 absorption at 21 cm (e.g., 
Garwood, Dickey 1989). There should also be a contribu- 
tion from the low-velocity dense clouds close to Sgr A. 
Some of them arc reported to be within about 50 pc of 
Sgr A by Whiteoak and Gardner (1978). 

5.1.4. +50 kmT x Wing 

A candidate for the H^" absorption wing at +50 km s _1 
is the "+50 km s _1 clouds," a complex of giant molec- 
ular clouds within about 10 pc of the Galactic nucleus 
(Giisten et al. 1981; Giisten, Henkel 1983). In the model 
of the large-scale structures in the central 10 pc of our 
Galaxy constructed with recent high spatial resolution 
radio observations (Coil, Ho 2000; Wright et al. 2001), 
the Sgr A East non-thermal radio source is impacting the 
"+50 km s _1 clouds" at the far side of the Galactic cen- 
ter. If the H^" absorption wing at +50 km s _1 is associated 
with these giant molecular clouds, it places GC IRS 3 be- 
yond Sgr A, which then jeopardizes the membership of 
GC IRS 3 in the central star cluster. Instead, to account 
for the absorption feature of CO in the near-infrared at 
4.7 fim, Geballe, Baas, and Wade (1989) proposed that the 
particular cloud occulting GC IRS 3 could be a part of the 
"circumnuclear disk" delineated by Giisten et al. (1987) 
in their HCN map. The "circumnuclear disk" is a com- 
pact (~ 3 pc) clumpy torus, a reservoir supposedly feeding 
mass to the dynamical center of our Galaxy, Sgr A* . The 
+50 km s" 1 component is barely seen in GCS 3-2, which 
suggests that the absorber of GC IRS 3 is a compact lo- 
cal structure at Sgr A, and thus argues in favor of the 
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"circumnuclear disk" origin. 

5.2. Special Characteristic of the Number Density 

Lines of sight toward the Galactic center sample both 
dense and diffuse clouds along the long pathlcngth. About 
30% of the visual extinction toward the Galactic center is 
believed to arise in dense clouds, and the rest in diffuse 
clouds (e.g. Whittet et al. 1997). For instance, the diffuse 
and dense cloud extinctions to GC IRS 3 estimated from 
the optical depths of hydrocarbon and water ice in the 
3 /jm region (Chiar et al. 2002) are Ay = 36 and Ay = 
11, respectively. 

The optical depth of Hg~, however, does not scale with 
the visible extinction. In contrast to most other molecules, 
the number density of Hjj" is constant and independent of 
the density of the cloud as long as the number density of 
H2 relative to that of the destroyer X of H^~, n(H2)/n(X), 
is constant (figure 6). A simple analysis using a steady 
state chemical kinetics yields a number density of 
?i(H^) = (C/fcco)[«-(H 2 )/ri(CO)] for a typical dense cloud 
where carbon atoms are mostly in the form of CO, and 
n(H^) = (C/fc c )[n(H 2 )/n(e)] for a typical diffuse cloud 
where carbon atoms arc in the form of C + (( is the cos- 
mic ray ionization rate, and fcco and k e are the rate con- 
stants of the ion-neutral reaction between and CO and 
of dissociative electron recombination of Hjj", respectively; 
McCall et al. 1998a,b). If typical values of C ~ 3 x 10" 17 
s -1 ) k C o ~ 2 x 10~ 9 cm 3 s _1 , fc c ~5x 10~ 7 cm 3 s _1 are 
used, and [n(H 2 )/n(CO)] and [n(H 2 )/n(e)] are both as- 
sumed to be 7 x 10 3 , we obtain n(ILj~) ~ 1 x 10~ 4 cm~ 3 
in dense clouds and ~ 4 x 10~ 7 cm~ 3 in diffuse clouds. 
These constancies of the number density make H^" nearly 
a direct indicator of the dimension of the absorbing clouds. 
In other words, regardless of how high or low the cloud 
density is, the Hj 4 " column density is simply proportional 
to the physical dimension of the cloud along the line of 
sight. This simple analysis well explains the observed Hjj" 
column densities in dense clouds (McCall et al. 1999). 

However, a problem appears when this model is applied 
to diffuse clouds because of uncertainties of the constants 
used in the calculation. Since the observed H^" column 
densities in dense and diffuse clouds are comparable, the 
large factor of 250 difference in n(H^") leads to the same 
large difference in the cloud dimension, which is difficult to 
accept (McCall ct al. 1998a,b; Geballe et al. 1999; McCall 
et al. 2002). Currently, at least one of the three assumed 
values, [C, fee, and n(e)] is under suspicion. For instance, 
measurements of k e vary by more than one order depend- 
ing on the experimental techniques (e.g. Larsson 2000). It 
is likely that the estimate of n(R^) in diffuse clouds will 
be increased by one order of magnitude or so when the 
uncertainty in these three parameters is resolved; n(ihj~) 
will then give a reasonable cloud dimension. 

5.3. in Diffuse Clouds 

The sharp and weak H2CO lines toward Sgr A at 
— 140 km s^ 1 , —60 km s _1 and km s _1 , and the weak- 
ness of the CO absorption at these velocities in GC IRS 3 
indicate that the clouds responsible for those absorptions 



are diffuse. Geballe, Baas, and Wade (1989) estimated the 
CO column density in the cloud at —60 km s" 1 velocity 
to be 1 x 10 17 cm -2 with a cloud temperature of 17 K. 
If we use [H 2 ]/[CO] = 7 x 10 3 , we obtain N(B 2 ) = 7 x 
10 20 cm" 2 , though this value might be better taken as a 
lower limit because of the possible saturation in the CO 
fundamental lines. Such a low hydrogen column density 
is consistent with diffuse clouds. 

The observed H3" column density in each of these ve- 
locity components of GC IRS 3 is in the range iV^Hg") 
= (2.4-5.6) x 10 14 cm" 2 . If the Iht number density in 
diffuse clouds given above (4 x 10 -7 cm -3 ) is assumed, 
we obtain unreasonably high cloud dimensions of 200 to 
500 pc. The enigma of H^" chemistry in the diffuse inter- 
stellar medium needs to be solved before we discuss the 
cloud dimension any further. 

The large pedestal component of GC IRS 3 is also in- 
terpreted to be due to diffuse clouds, because of the ab- 
sence of strong CO or H2CO absorption. A clarification of 
the nature of this cloud is an interesting future problem. 
There exists a possibility that the decomposition of the 
observed features into velocity components is not unique 
in view of the significant noise in the spectrum and the 
limited spectral resolution. 

5.4- H3 in Dense Clouds 

The strong and saturated CO absorption toward 
GC IRS 3 and H 2 CO for Sgr A at the velocity of 
+50 km s _1 wing clearly represents dense clouds (figure 
5). It is remarkable that the ITj~ absorption lines starting 
from the J = 1 levels are weak at this velocity (figure 2). 
This must be due to a small pathlcngth of a very dense 
cloud in front of the infrared source. From the equiv- 
alent width of the R(l, I) 1 line integrated over +20 to 
+60 km _1 , N(R+) J=1 = 2.4 xlO 14 cm" 1 is obtained; 
this corresponds to a dense cloud pathlength of 1.1 pc if 
n(Hj) ~lx 10 -4 cm" 3 estimated in subsection 5.2 is as- 
sumed. The concentration of the H^" in the (3, 3) level 
at the velocity would increase the total column density 
and the cloud dimension; however, we leave this popula- 
tion not included because we do not have sufficient knowl- 
edge about either the nature of the cloud that holds Hjj" 
in the (3, 3) level or the mechanism that populates the 
metastable states, which we discuss separately in the next 
subsection. 

The cloud dimension of the +50 km s~ J wing might not 
be incompatible with the "+50 km s _1 clouds" interpre- 
tation in terms of the physical parameters presented by 
Glisten and Hcnkcl (1983). They estimate the hydrogen 
column density of M-0.02-0.07, the closest "+50 km s _1 
cloud" to GC IRS 3, to be 2V(H 2 ) = 7 x 10 24 cm~ 2 , which 
gives n(H 2 ) = 2 x 10 6 cm -3 if a cloud dimension of 1 pc 
is assumed. The angular diameter of M-0.02-0.07 is 2' in 
Giisten and Henkel (1983), which is 5 pc at the assumed 
8 kpc distance. The apparent smaller pathlcngth mea- 
sured by could be reconciled if the line of sight toward 
GC IRS 3 is off-centered from the cloud core. 

On the other hand the physical model of the "cir- 
cumnuclear disk" proposed by Marshall, Lasenby, and 
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Harris (1995) gives the radial thickness of the annular disk 
to be 0.5 pc with a hydrogen density of n(H.2) — 10 6 cm~ 3 . 
The factor 2-3 difference could be filled by the geometry 
of GC IRS 3 and the disk. 

5.5. in Metastable State 

The tentative detection of the i?(3, 3)' transition at 
3.534 /im for the first time in interstellar space would be 
a breakthrough if confirmed by more observations; it in- 
troduces a new dimension in the studies of interstellar . 
The detection of i?(3, 3)' and non-detection of other tran- 
sitions starting from J > 1 levels, such as R(2, 1)", i?(2, 
2)', Q(2, l) 1 , Q(3, 0) etc., is reasonable since only the 
(3, 3) rotational level is metastable; spontaneous emission 
from this level to lower levels is forbidden by the ortho- 
para selection rule and the absence of the (2, 0) rotational 
level by the Pauli principle, as shown in figure 7 (Pan, Oka 
1986). On the other hand, H+ in the (2, 1), (2, 2), (3, 0), 
(3, 1), and (3, 2) levels decays to lower levels with lifetimes 
of 20.4 d, 27.2 d, 3.8 hr, 7.9 hr and 15.8 hr, respectively 
(Neale et al. 1996) through centrifugal distortion-induced 
rotational transitions initially proposed for NH3 (Oka et 
al. 1971). Since detailed ab initio quantum chemical cal- 
culations of H^" have been extensively carried out, these 
lifetimes arc accurate and reliable. All ortho-H^j" which 
are chemically produced or collisionally pumped into the 
(4, 3), (3, 0), and (3, 3) levels accumulate in the (3, 3) level 
until it is collisionally deexcited to lower levels. Similarly, 
ortho-H^ produced in (5, 3), (5, 0), (6, 6), (6, 3), (7, 6), 
(7, 3), (7, 0) etc. accumulate in the (6, 6) metastable 
level. Metastable levels are also expected for para-H^, for 
which (4, 4) and (5, 5) may be excessively populated. H3 
in the (4, 4) level can spontaneously decay to the (3, 1) 
level, but its lifetime is long (11 yr). 

Hjj" in metastable levels relaxes to lower levels by col- 
lisions with H 2 . Unlike NH 3 in which ortho «-> para 
transitions are forbidden for collision-induced transitions 
(Cheung et al. 1969), may relax to lower levels of dif- 
ferent spin modification. This is because a collision of 
with H2 is a chemical reaction and the scrambling of 
protons may change ortho-H^ into para-H^ or vice versa, 
although some selection rules for nuclear spin still remain 
(Uy et al. 1997; Cordonnicr et al. 2000). 

There are two ways that the (3, 3) level could be sig- 
nificantly populated compared with the lowest J = 1 lev- 
els. First, it may be simply because of high temperature. 
The (3, 3) level is higher than the lowest (1, 1) level by 
251.23 cm -1 (361 K) and the populations of the two levels 
become equal at T = 234 K. Populations in excited non- 
metastable levels may still be low even at high temper- 
ature because of fast spontaneous emission, if the cloud 
density is lower than the critical densities of the levels. 
The critical densities are on the order of 10 4 cm~ 3 for (2, 
1) and (2, 2) and 10 6 cm" 3 for (3, 0), (3, 1) and (3, 2). 
Thus our observation of the i?(3, 3)' transition and non- 
detection of other transitions starting from J > 1 levels 
indicate low density. If the metastable state is positively 
confirmed to be thermally populated, then the i?(3, 3)' 
transition and other absorption lines from metastable lev- 
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els will make an excellent probe to efficiently isolate the 
warm interstellar gas. These absorption lines are unique in 
that they are sensitive exclusively to the high temperature 
clouds (T ~ 200 K) along the line of sight with minimum 
contamination by other cold clouds. 

Second, the excess population at (3, 3) could occur at 
low temperature with a subtle balance of the rate of col- 
lisional deexcitation and destruction of H3". If the for- 
mer rate, fcH?i(H;j~)n(H2), is much faster than the latter, 
fcxi(H| )n(X), at low temperature will all relax to the 
lowest J = 1 levels during its lifetime. The excess pop- 
ulation in (3, 3) could occur only when those rates are 
comparable, that is, /ch^H^) ~ fcx"(Hx)- This cannot 
happen in dense clouds where CO is the main destroyer 
of , since fcn ~ fcco an d n(H 2 ) 3> n(CO). For diffuse 
clouds this may be possible if n(H2)/n(e) is not much 
larger than k e /kn ~ 250. For this to happen a high elec- 
tron density and high k c are preferred. It is interesting 
to note that the latter requirement is opposite to what 
is needed to reconcile the enigma of the observed high 
column densities of rpj~ in a diffuse interstellar medium. 

Although our detections of the i?(3, 3) z line look fairly 
convincing for both GC IRS 3 and GCS 3-2, we note that 
both lines are broad and without a sharp velocity compo- 
nent. In addition, the spectral line toward GC IRS 3 is 
enigmatic in that the five spectral lines starting from the 
J = 1 level are not strongly observed at the velocity of the 
i?(3, 3) z line. A simple model calculation shows that it is 
hard to populate the (3, 3) level alone without populating 
the (1, 0) and (1, 1) levels. It is intriguing that the line 
profile of i?(3, 3)' is similar to those of the strongly sat- 
urated CO line both in its wavelength and widths (figure 
5). Such CO lines indicate dense clouds where will 
be collisionally cooled rapidly to lower levels. We may 
be sampling low-density warm gas with a long pathlcngth 
surrounding such dense cloud. 

The high population in the (3, 3) level makes the prob- 
lem of the unexpectedly high abundance in diffuse clouds 
even harder to reconcile, since other metastable states 
may also be significantly populated. The need of a high 
value of k e for the second mechanism to populate (3, 3) 
is also an interesting twist related to the enigma. We 
will attempt to observe more objects and more spectral 
lines starting from metastable levels, and to carry out a 
detailed analysis of this problem. 

6. Summary 

We have presented an absorption-line survey of 
toward two Galactic center sources: GCS 3-2 in the 
Quintuplet cluster and GC IRS 3 near Sgr A*. Six 
lines were detected for each source, of which three were 
newly detected. In particular, the absorption of H;j~ orig- 
inating from the (J, K) = (3, 3) metastable state was 
tentatively detected for the first time in interstellar space. 
The observed absorption lines show intriguing line 
profiles, indicating at least four velocity components of 
clouds along the line of sight. The velocities well match 
those of H 1, CO, and H 2 CO reported earlier in the radio 
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and infrared. The Hg velocity components at 0, —60, and 
— 110 - —140 km s -1 , which well match those of sharp 
and weak CO as well as with H2CO, are inferred to be in 
the diffuse interstellar medium of intervening spiral arms, 
while the +50 km s _1 component in GC IRS 3 should be 
associated with the local structure of the Galactic nucleus. 
Surprisingly, not much H^" is found at +50 km s _1 where 
very strong and saturated CO absorption is observed in 
GC IRS 3. This shows the marked contrast between CO 
and as astrophysical probes; the number density of 
CO is proportional to the cloud density, while that of H^ 
is independent of the cloud density. 

The tentative detection of the i?(3, 3)' line provides 
observational evidence for the metastability of the (3, 3) 
level, which had been theoretically expected. It suggests 
that other metastable J = K levels may also be populated. 
While the observed H^ spectra reveal the great richness of 
the sight lines toward the Galactic center, there are many 
loose ends in our interpretation of the spectra. Further 
studies on the H^ spectrum as well as on CO and H2 
spectroscopy toward the Galactic center will be attempted 
to further clarify this situation. 
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Table 1. Major absorption lines of HJ" vi <— targeted in 
the survey. 



Transition 


ALab 




Coverage 


W x 


Mevel 




[Mm]* 


[D 2 ]t 




[10 


5 /im] 


[10 15 


cm 2 ] 










GCS 3-2 


GC IRS3 


GCS 3-2 


GC IRS3 


i?(3, 3)" ... 


3.426974 


0.0071 


No 










i?(3, 3) ; ... 


3.533666 


0.0191 


Yes 


2.8 


3.2 


1.0 


1.1 


R(2, \) u . . . 


3.538424 


0.0182 


Yes 










R(2, 2)" ... 


3.542158 


0.0094 


Yes 










r(2, iy ... 


3.615924 


0.0044 


No 










R(2, 2)' ... 


3.620473 


0.0177 


No 










R(l, 1)" ... 


3.668083 


0.0158 


Yes 


lit 


8.2* 


2.5 § 


1.9 § 


J2(l, 0) ... 


3.668516 


0.0259 


Yes 


11* 


8.2* 


1.2$ 


0.94§ 


iy ... 


3.715479 


0.0141 


Yes 


5.111 


4.811 


2.3 


2.2 


0(3, 3) ... 


3.903967 


0.0065 


No 










0(2, 2) ... 


3.914406 


0.0086 


Yes 










s \ / n 1 Nil 

0( 2 , 1) ■ ■ 


3.916979 


0.0043 


"XT' 

Yes 










0(1, i) ... 


3.928625 


0.0128 


Yes 


4.5 


5.1 


2.1 


2.5 


0(1, o) ... 


3.953000 


0.0254 


Yes 


5.2 


5.0 


1.3 


1.2 


0(2, 1)' ... 


3.971073 


0.0167 


Yes 










0(3, 0) ... 


3.985528 


0.0249 


Yes 










0(3, 1)" .. 


3.987028 


0.0217 


Yes 










P(l, 1)... 


4.069524 


0.0086 


No 










iVtotal [10 i5 


cm"' 1 ] 










4.6 


4.4 


Note. The 


laboratory data 


are from 


McKellar 


and 









Watson (1998) and Lindsay and McCall (2001). The 
velocity-integrated equivalent widths and the column den- 
sities are summarized for the detected lines. 
* Lindsay and McCall (2001). 
t J. K. G. Watson, private communication, 
t Total equivalent width of R(l, l) u and R(l, 0). 
§ Column density of the each level is derived on the 
assumption that the populations in the ortho- (1,0) and 
para- (1, 1) levels are divided into about 2 : 1, which is 
implied by R(l, l) 1 , 0(1, 1) and 0(1, 0) transitions. 
II Equivalent width of R(l, l) 1 is calculated from the 
Gaussian decomposition parameters. For other tran- 
sitions Wx are obtained by directly integrating the 
absorption features. 
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Name 


R.A. 


Dec. 


I 


b 


L 


Exposure 


Grating* 


Standard 






(J2000) 


(J2000) 


[°] 


[°] 






ECH 


XDP 


Name 


Spe. 


GCS 3-2T. . . 


17:46:14.8 - 


-28:49:41 


+0.16 - 


-0.06 


2.7 


120 


8350 


6100 


HR 7528 


B9.5 IV 














240 


4400 


5200 


HR 7528 


B9.5 IV 


GC IRS 3*... 


17:45:39.9 - 


-29:00:24 


-0.06 - 


-0.04 


5.3 


360 


8350 


6100 


HR 7924 (a Cyg) 


A2 lac 














1440 


4400 


5200 


HR 7924 (a Cyg) 


A2 Iae 


: "ECH" and 


"XDP" denote 


the angle 


of cchellc 


and 















cross-dispersing gratings in the instrumental unit. 

t Coordinate and magnitude are from Nagata et 

al. (1990). 

* Coordinate and magnitude are from Blum, Scllgren, 
and DePoy (1996). 
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Table 3. Velocity-resolved components of R(l, l) 1 , Q(l, 1) 
and Q(l, 0). 



NcLIIlC 




FWHM 




W\ 

' " A 






1 v level 




1 v total 




[km s -1 ] 


[km s -1 ] 




flO -5 uml 






no 14 cm~ 2 l 




[10 14 cm -2 ]* 








K{L, 1J 




U) 










GCS 3-2. . . 


-165 


15 


0.15 






0.7 






1.3 




-107 


50 


1.8 


2.0 


1.2 


8.3 


9.4 


2.9 


16 




-51 


38 


1.9 


1.7 


2.1 


8.6 


7.9 


5.1 


17 




-2 


31 


1.2 


1.5 


1.7 


5.6 


7.0 


4.1 


11 


GC IRS 3. . . 


-140 


22 


0.64 




0.30 


2.9 




0.51 


5.6 




-56 


15 


0.27 




0.46 


1.2 




1.1 


2.4 




-23 t 


116 


3.7 




3.0 


17 




7.2 


32 




-5 


12 


0.50 




0.73 


2.3 




1.7 


4.3 



* Based on the velocity independent conversion factor of 
R(l, I) 1 to total column density derived from table 1. 
t Pedestal component. 
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Fig. 1. Wavelength coverage of the IRCS echelle mode in the 
3 (im region. The entire L-band (2.84-4.18 (im) is covered by 
six standard settings (LA - , LAo, LA" 1 ", LB~, LBo, LB + ; 
shown in green and blue rectangles), but 70% is covered by 
three settings. "ECH" and "XDP" denote the angles of the 
echelle and the cross-disperser gratings. The grating angle 
can be seamlessly modified so as to be best suited to indi- 
vidual projects. The two settings which we employed for the 
observation are shown by the bold-line enclosures. The posi- 
tions of the major Hg" lines are indicated by sticks with the 
heights representing the relative intensities calculated with 
T = 300 K. The atmospheric transmission curve is overlaid 
to illustrate the region where the interference of the telluric 
atmospheric absorption is severe. 



Fig. 6. Schematic of the number density variation of Hg 
with respect to hydrogen along with that of major molecules 
and ions in the interstellar medium. The number density 
of Hg" is independent of n(H), in contrast with the other 
molecules. We define "dense" and "diffuse" clouds by the 
form of carbon atoms. The form of carbon is critical, since 
it defines the destruction mechanism of Hj, and hence the 
number density of the molecular ions. 



Fig. 2. Total detected Hg absorption spectra toward 
GCS 3-2 (left column) and GC IRS 3 (right column). 
Kinematic components arc indicated in the LSR velocity by 
dotted lines (See figure 4). Note that the 3.668 fim line profile 
appears to be different from the others because it consists of 
a close doublet of R(l, 0) and R(l, l) u separated by Av = 
35 km s — 1 . The feature is lined up with other absorption lines 
at R(l, 0) 3.6685 (im. The model atmospheric transmission 
curves are shown to discriminate the genuine detections from 
the possible artifacts of the poor cancellation of the telluric 
absorption lines. 



Fig. 3. Same with figure 2, but for negative detections. 
The periodic patterns seen at 3.97-3.99 (im are instrumen- 
tal fringes left unprocessed. Note that the dips at 3.534 (im 
in the bottom panels arc i?(3, 3)' positively detected. 



Fig. 4. Gaussian decomposition of the kinematic compo- 
nents for GCS 3-2 (top) and GC IRS 3 (bottom). At least 
four components are found in ortho- Q(l, 0) and para- R(l, 
l) 1 and Q(l, 1) H^~ absorption lines at common velocities. 
The residual of the fitting is shown at the bottom of each 
panel. The decomposition parameters are summarized in ta- 
ble 3. 



Fig. 5. Comparison of the line profiles of Hg" at R(l, 1)' 
and R(3, 3) along with those of other species toward the 
same objects. GCS 3-2 is on the left, and GC IRS 3 on the 
right. The absorption spectra from previous studies are shown 
in arbitrary units. The dotted lines represent the velocity 
components of the R(l, l) 1 transition in the previous figure. 
Three velocity components in GC IRS 3 at —140 km s — , 
-60 km s" 1 , and km s -1 match well with CO R(2), H 2 CO 
and H I absorption, indicating the same absorbing clouds. 
The deep absorptions of CO and H2CO at +50 km s -1 to- 
ward the Galactic center do not correspond to strong com- 
ponents in the H^~ spectrum, which might be accounted for 
by the compact structure of the absorbing cloud. Note the 
considerable difference in the line shape of R(3, 3)' from R(l, 
l) 1 for both GCS 3-2 and GC IRS 3. The line profile of the 
i?(3, 3)' in GC IRS 3 is more or less similar to that of CO 
absorption at 4.7 fira. 



Fig. 7. Rotational energy levels of Hg in the ground state 
and in the V2 vibrationally excited state. Hg*" has two vibra- 
tional modes of which only V2 is infrared active, and targeted 
in our observation. G = \k — Z2I is used instead of K in the 
diagram, where I2 is a vibrational angular momentum of 1^2- 
G is a better quantum number than K where I2 7^ (see 
McCall 2001 for details). The letters u and I marked in the 
1/2 state denote the levels of the same G number, but with 
different combinations of k and Z2. The major transitions of 
V2 <— that appear in the 3 to 4 fim regions are indicated 
by the vertical connecting arrows with labels. Note that for 
the ground vibrational state, the J = 2n and G = levels are 
forbidden by the Pauli principle (shown in dotted bars). The 
lowest levels of ortho- (J, G) = (1, 0) and para- (J, G) = (1, 
1) Hg" are marked with thick solid bars. The selection rules 
allow only radiative relaxation between two levels that satisfy 
A J = 0, ±1 and AG = ±3. The spontaneous transitions are 
indicated with connecting lines of energy levels in the ground 
state. The transition between G = 2 <-> 1 can be understood 
to be G = ±2 <-> +1. The "+" and "— " signs in the ground 
state denote the parity of the energy levels. The parity rule, 
+ «-> — , is automatically satisfied by AG = ±3. The (J, G) 
= (3, 3) and (5, 5) are disconnected from any lower levels 
by the selection rule, making the two levels metastable (Pan, 
Oka 1986). The (4, 4) level can relax to (3, 1), but has a long 
lifetime. 
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v LSR [km/s] 
-500 -250 250 500 



= 0.95 



0.9 
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